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Abstract — We propose the QoS-aware BS-selection and the cor- 
responding resource-allocation schemes for downlink multi-user 
transmissions over the distributed multiple-input-multiple-output 
(MIMO) links, where multiple location-independent base-stations 
(BS), controlled by a central server, cooperatively transmit data to 
multiple mobile users. Our proposed schemes aim at minimizing 
the BS usages and reducing the interfering range of the distributed 
MIMO transmissions, while satisfying diverse statistical delay-QoS 
requirements for all users, which are characterized by the delay- 
bound violation probability and the effective capacity technique. 
Specifically, we propose two BS-usage minimization frameworks to 
develop the QoS-aware BS-selection schemes and the correspond- 
ing wireless resource-allocation algorithms across multiple mobile 
users. The first framework applies the joint block-diagonalization 
(BD) and probabilistic transmission (PT) to implement multiple 
access over multiple mobile users, while the second one employs 
time-division multiple access (TDMA) approach to control multi- 
ple users' links. We then derive the optimal BS-selection schemes 
for these two frameworks, respectively. In addition, we further 
discuss the PT-only based BS-selection scheme. Also conducted is 
a set of simulation evaluations to comparatively study the average 
BS-usage and interfering range of our proposed schemes and to 
analyze the impact of QoS constraints on the BS selections for 
distributed MIMO transmissions. 

Index Terms — Distributed MIMO, broadband wireless net- 
works, statistical QoS provisioning, wireless fading channels. 

I. Introduction 

TO increase the coverage of broadband wireless networks, 
distributed multiple-input-multiple-output (MIMO) tech- 
niques, where multiple location-independent base stations (BS) 
cooperatively transmit data to mobile users, have attracted 
more and more research attentions |Q]-|l3]. In particular, the 
distributed MIMO techniques can effectively organize multiple 
location-independent BS's to form the distributed MIMO links 
connecting with mobile users Like the conventional centralized 
MIMO system (H-||6l, the distributed MIMO system can 
significantly enhance the capability of the broadband wireless 
networks in terms of the quality-of-service (QoS) provision- 
ing as compared to the single antenna system. However, the 
distributed nature for cooperative multi-BS transmissions also 
imposes many new challenges in wide-band wireless com- 
munications, which are not encountered in the centralized 
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MIMO systems. First, the cooperative distributed transmissions 
cause the severe difficulty for synchronization among multiple 
location-independent BS transmitters. Second, as the number 
of cooperative BS's increases, the computational complexity 
for MIMO signal processing and coding also grow rapidly. 
Third, because the coordinated BS's are located at different 
geographical positions, the cooperative communications in fact 
enlarge the interfering areas for the used spectrum, thus dras- 
tically degrading the frequency-reuse efficiency in the spatial 
domain. Finally, many wide-band transmissions are sensitive to 
the delay, and thus we need to design QoS-aware distributed 
MIMO techniques, such that the scarce wireless resources can 
be more efficiently utilized. 

Towards the above issues, many research works on dis- 
tributed MIMO transmissions have been proposed recently. 
The feasibility of transmit beamforming with efficient syn- 
chronization techniques over distributed MIMO link has been 
demonstrated through experimental tests [2], suggesting that 
complicated MIMO signal processing techniques are promising 
to implement in realistic systems. For the centralized MIMO 
system, the antenna selection 0, is an effective technique 
to reduce the complexity, which clearly can be also extended 
to distributed MIMO systems for the BS selection. It can 
be expected that the BS-selection techniques can significantly 
decrease the processing complexity, while still achieving high 
throughput gain over the single BS transmission. Also, it is 
desirable to minimize the number of selected BS's through BS- 
selection techniques, which can effectively decrease the inter- 
fering range and thus improve the frequency-reuse efficiency 
of the entire wireless network. Most previous research works 
for BS/antennas selections mainly focused on the scenarios 
of selecting a subset of BS's/antennas with the fixed cardi- 
nality 0, 0, 0. However, it is evident that based on the 
wireless-channel status, BS-subset selections with dynamically 
adjusted cardinality can further decrease the BS usage. More 
importantly, how to efficiently support diverse delay-QoS re- 
quirements through BS-selection in distributed MIMO systems 
sill remains a widely cited open problem. 

To overcome the aforementioned problems, we propose the 
QoS-aware BS-selection schemes for the distributed wireless 
MIMO links, which aim at minimizing the BS usages and 
reducing the interfering range, while satisfying diverse sta- 
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tistical delay-QoS constraints. In particular, we develop two 
BS-usage minimization frameworks for distributed multi-suer 
MIMO transmissions. The first framework uses the joint block- 
diagonalization (BD) and probabilistic transmission (PT) for 
multiple access of multi-user over distributed MIMO links, 
while the second framework employs time-division multiple 
access (TDMA) techniques. We derive the optimal QoS- 
aware BS-selection and the corresponding resource allocation 
schemes for these two frameworks, respectively. We also 
discuss the PT-only based BS-selection scheme. Simulations 
are conducted for comparative analyses among the above BS- 
selection schemes. 

The rest of this paper is organized as follows. Section HI] de- 
scribes the system model for distributed MIMO transmissions. 
Section [HI] introduces the statistical QoS guarantees and the 
concept of effective capacity. Section HVl develops the joint BD 
and PT (BD-PT) optimization frameworks for QoS-aware BS- 
selections over multi-user distributed MIMO links and derives 
the corresponding optimal solution. Section IVl derives TDMA- 
based QoS-aware BS-selection scheme. Section [VI] conducts 
simulations to perform comparative analyses for our proposed 
schemes. The paper concludes with Section IVIII 

Notations: The operator | • | used on a real or complex 
number generates the absolute value; the operator | • | used for 
a set represents the cardinality of this set. We use boldface 
to denote matrices and vectors. For an X x Y matrix A, 
we denote by A(i,j) the element on the ith row and jth 
column; |A||_f denotes the Frobenius norm of A, where 

II A||f = JEliEU \Mhj)\ 2 - The operators and (•)+ 
generate the transpose and conjugate transpose, respectively. 
The operator 1/.) is the indication function. If the statement in 
the subscript is true, we have — 1; otherwise, 1/.) = 0. 

II. System Model 
A. System Architecture 

We concentrate on the wireless distributed MIMO system 
for downlink transmissions depicted in Fig.Q] which consists of 
Kfr s distributed BS's, K mu mobile users, and one central server. 
The mth BS has M m transmit antennas for m = 1, 2, . . . , K^, s 
and the nth mobile user has N n receive antennas for n = 
1,2,..., K mu . All distributed BS's are connected to the central 
server through high-speed optical connections. The data to be 
delivered to the nth mobile user, n = 1, 2, . . . , K mu , arrives at 
the central server with a constant rate, which is denoted by C n . 
Then, the central server dynamically controls these distributed 
BS's to cooperatively transmit data to the corresponding mobile 
users under the specified delay-QoS requirements. 

For multi-user downlink transmissions, the distributed BS's 
and the mobile users form the broadcast MIMO link for data 
transmissions. The wireless fading channels between the mth 
BS and the nth mobile user is modeled by an N n x M m 
matrix H„ m , where H n ,m(hj) is the complex channel gain 
between the ith receive antenna of nth mobile user and the 
jth transmit antenna of the mth BS. All elements of H JljTO 
are independent and circularly symmetric complex Gaussian 
random variables with zero mean and the variance equal to 
h n ,m, implying that H has continuous cumulative distribution 




Fig. 1. System model for wireless downlink distributed MIMO transmissions. 

function (CDF). Also, the instantaneous aggregate power gain 
of the MIMO link between the ?ith mobile user and the mth 
BS, denoted by 7„, m , is defined by 

^In.ra — ~7J ||H nim ||^ (1) 

Since the Frobenius norm of the channel matrix can effectively 
characterize the channel quality in terms of achieving high 
throughput [5], the aggregate power gain given in Eq. ([T} will 
play an important role in our BS selection design. We further 
define H„ = [H„,i H rlj 2 ■ ■ ■ H„,K bB ] as the CSI for the nth 
mobile user for n = 1,2,..., K mu . The matrix H„ follows the 
independent block-fading model, where H„ does not change 
within a time period with the fixed length T, called a time 
frame, but varies independently from one frame to the other 
frame. Furthermore, we define H = [H[ H2 • ■ • H^ b ] T , 
representing a fading state of the entire distributed MIMO 
system. 

In order to decrease the complexity and suppress the inter- 
fering range of the distributed MIMO transmission, the cental 
server dynamically selects a subset of BS's to construct the 
distributed MIMO link. Then, our design target is to minimize 
the average number of needed BS's subject to the specified QoS 
constraints. We suppose that each mobile user can perfectly 
estimate its CSI at the beginning of every time frame and 
reliably feed CSI back to the central server through dedicated 
control channels. Based on CSI H and QoS requirements, the 
central server then adaptively selects the subset of BS's and 
organizes them to transmit data to mobile users through the 
distributed MIMO links. 



B. The Delay QoS Requirements 

The central data server maintains a queue for the incoming 
traffic to each mobile user. We mainly focus on the queueing 
delay in this paper because the wireless channel is the major 
bottleneck for high-rate wireless transmissions. Since it is 
usually unrealistic to guarantee the hard delay bound over the 
highly time-varying wireless channels, we employ the statis- 
tical metric, namely, the delay-bound violation probability, to 
characterize the diverse delay QoS requirements. Specifically, 
for the nth mobile user, the probability of violating a specified 

(n) 

delay bound, denoted by Dl, , cannot exceed a given threshold 
£„. That is, the inequality 

Pr{£>„ >D$} <f», n = l,2,...,N mu , (2) 
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needs to hold, where D n denotes the queueing delay in the nth 
mobile user's queueing system. 

C. Performance Metrics and Design Objective 

We denote by L the cardinality of the selected BS subset 
(the number of selected BS's) for the distributed MIMO trans- 
mission in a fading state. Then, we denote the expectation of 
L by L and call it the average BS usage. As mentioned in 
Section Hi- Al our major objective is to minimize L through dy- 
namic BS selection while guaranteeing the delay QoS constraint 
specified by Eq. (ffjl. We will also evaluate the average interfer- 
ing range affected by the distributed MIMO transmission. The 
instantaneous interfering range, denoted by A, is defined as the 
area of the region where the average received power under the 
current MIMO transmission is larger than then certain threshold 
denoted by er| h . The average interfering area is then defined as 
the expectation Eh {^4} over all H. Clearly, minimizing L can 
not only reduce implementation complexity, but also decrease 
the average interfering range affected by the transmit power. 

D. The Power Control Strategy 

The transmit power of our distributed MIMO system varies 
with the number of selected BS's. In particular, given the 
number L of selected BS's, the total instantaneous transmitted 
power used for distributed MIMO transmissions is set as a 
constant equal to Vl- Furthermore, Vl linearly increases with 
L by using the strategy as follows: 

V L = V iei + K(L-l), L=l,2,...,K hs , (3) 

where V T el > is called the reference power and k > 
describes the power increasing rate against L. Also, we define 
Vl = for L = 0. The above power adaptation strategy 
is simple to implement, while the average transmit power 
can be effectively decreased through minimizing the average 
number of used BS's. In addition, Eq. (0 can upper-bound the 
instantaneous interferences and the interfering range over the 
entire network. 

III. Effective Capacity Approach for Statistical 
Delay-QoS Guarantees 

In this paper, we apply the effective capacity approach JS], 
|9l , ifTTI . ifTTl to integrate the constraint on delay-bound 
violation probability given by Eq. (0 into our BS selec- 
tion design. Consider a stable discrete-time queueing system 
with the stationary time-varying arrival-rate and departure-rate 
(service-rate) processes. The asymptotic analyses based on the 
large deviation principal Q, JH) show that under the sufficient 
conditions, the probability that the queue-length, denoted by 
Q, exceeding a given bound Q t h can be approximated by 

Pr{Q > Qth} « e- eQ <\ (4) 

where 9 > is a constant called QoS exponent. It is clear that 
the larger (smaller) 9 implies the lower (higher) queue-length- 
bound violation probability. 

By using 9, the delay-bound violation probability can be 
approximated (8) by 

Pr{Z?> Ah} « e- eZ?Dt >\ (5) 



for the constant rate C. When the arrival rate is not time- 
varying, the approximation in Eq. needs to replace C with 
effective bandwidth [7], |8| function of the arrival rate process, 
which is defined as the minimum constant service rate required 
to guarantee QoS exponent 9. 

Then, to upper-bound Pi{D > Ah} with a threshold £, 
using Eq. 0, we get the minimum required QoS exponent 9 
as follows: 

e=-m. (6) 

CAh 

Consider a discrete-time arrival process with constant rate C 
and a discrete-time time-varying stationary departure process, 
denoted by R\k\, where k is the time index. In order to 
guarantee the desired 9 determined by Eq. (|6), the statistical 
QoS theory [7|, [8| shows that the effective capacity C(9) of 
the service-rate process R[k] needs to satisfy 

C(9) = C, (7) 

given the QoS exponent 9. The effective capacity function is 
defined in [8] as the maximum constant arrival rate which can 
be supported by the service rate to guarantee the specified 
QoS exponent 9. If the service-rate sequence R[k] is time 
uncorrected, the effective capacity can be written IfTTl as 

C(9)^- l -\o g (^{e- dR ^}) 1 (8) 

where E{-} denotes the expectation. 

In our distributed MIMO system, the BS selection result 
is designed as the function determined by the current CSI. 
Thus, the corresponding transmission rate (service rate) is time 
independent under the independent block-fading model (see 
Section Ill-Al l. Then, applying Eqs. ©-([7]), the delay QoS 
constraints given by Eq. (fJJ can be equivalently converted to: 

E H {e- e " fl " _ e -0nC„j < 0; )n= i ) 2,...,JV ora , (9) 

where 9 n = — log(£„)/ / (C , n A^h' ) ) and Eh{-} denotes the 
expectation over all H. 

IV. Joint Block-Diagonalization and Probabilistic 
Transmission Based Base-Station Selection 

As discussed in Section [Tl] based on CSI H and QoS 
requirements, the central server will adaptively select the 
subset of BS's and organizes them to transmit data to mo- 
bile users through the distributed MIMO links. Given the 
cardinality L of the desired BS subset in a fading state, we 
denote by £Il the set of indices of selected BS's, where 
Ml = {iL,i,iL,2, ■ ■ ■ and i L ,e e {1, 2, . . . , K hs } for 

I = 1,2, ... ,L. Note that once a BS is selected, we use all 
its transmit antennas for data transmissions. For the specified 
L, we use Ul = {njj.i, n u,2, ■ ■ ■ , nu,u} to denote the set of 
active users, picked by the central server, which can receive the 
data in this fading state, where U is the cardinality of Ul- For 
presentation convenience, we use 9JIl = (£Il,Ml) to represent 
a specific transmission mode (or mode in short). Moreover, we 
term the pairs with U > 2 for Ul as multi-user modes, and 
term the pairs with U = 1 as single-user modes. 
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Given L, Ql and Ul, the channel matrix of the nth mobile 



user for n G Ul, modeled by H[™\ is determined by 



< ^ [H 



H„ 



(10) 



Mi; 



matrix. Furthermore, 

(n) 



where is an iV„ x (£^ efi 

we use Yq l to denote the power gain matrix for under 
the given f^, where 



r£'(i,j) = E H i H^(i,j) 



fixing f^i 



(11) 



The physical-layer signal transmissions can be modeled by 



(n) 



_(n) 



n G Wl, 



where s^j represents the ith user's input signal vector for the 
MIMO channel , y$ is the signal vector received by 
the nth user, and is the complex additive white Gaussian 
noise (AWGN) vector with unit power for each element of this 
vector. In this section, we employ the block-diagonalization 
(BD) technique l[T4ll to implement multiple access for multi- 
user modes in our QoS-aware BS-selection framework. 

For dynamic BS selections in distributed MIMO trans- 
missions, L and 50? l are both functions of CSI and QoS 
requirements. Then, we need to answer the following questions: 
(i) Given L, which transmission mode will be used for single- 
user and multi-user modes, respectively? (ii) When do we use 
single-user or multi-user modes? (iii) For a specific multi- 
user mode, how do we quantitatively allocate the wireless 
resources across multiple mobile users under the BD based 
transmissions? (iv) Which L will be selected for distributed 
MIMO transmissions in each fading state to decrease the 
average BS-usage while satisfying the QoS requirements? 

Clearly, we can not examine all combinations of {Q,l,Ul) to 
minimize the BS usage due to the too high complexity. Then, in 
Section ITY-A1 we develop the heuristic algorithms to efficiently 
select VJIl for the specified L in multi-user transmission 
modes. In Section IIV-B1 we determine how to select DJIl in 
single-user transmission mode. Based on schemes developed 
in Sections IIV-AI and IIV-B1 we further answer questions (iii) 
and (iv) through formulating and solving the joint BD-PT based 
BS-usage minimization problem in Sections IIV-CI and IIV-DI 



A. Selection of^fflL in Multi-User Transmission Modes 

In each fading state, we pick K^ s multi-user transmis- 
sion modes as candidates for distributed MIMO transmis- 
sions. These K^ s transmission modes corresponds to L = 
1,2,..., Kha, respectively, representing different levels of BS 
usages. As mentioned previously, the derivation of global 
optimal selection strategy in terms of minimizing the average 
BS usage is intractable, since the complexity of examining all 
possible DJXl = (SIlMl) is too high. Therefore, for a given L, 
we determine 9JTl through a two-step method. We first propose 
the priority BS-selection to determine the BS subset £1^. Then, 
based on the selected fli, we derive Ul through a joint channel- 
priority user-selection process. 

A. 1. Priority BS-Selection to Determine CIl 



01. Let * : = {1, 2, . . . , K hs }, * := 0, and t = |*[; ! Initialization 

02. j : = 1. ! Start selection with User 

03. WHILE (I < L) ! Iterative selections until L BS's are selected 

04. m* = argmin mSS {7^ (j)im }. 

! User ir(j) selects the BS with the largest aggregate power gain. 

05. * := * U {m*}, * := * - {m*}, and £ := £ + 1. 

! Update *, and £. 

06. IF j = K mu , then j := 1; ELSE j := j + 1. 

! Let next user with lower priority to select BS. 

07. END _ 

08. Q L := ¥. ! Complete the BS selection and get Q L . 

Fig. 2. The pseudo codes to determine Ql m eacn fading state by using the 
priority BS-selection algorithm for the multi-user transmissions. 

Consider any fading state H. The nth user's global maximum 
achievable transmission rate is attained when all BS's are used 
and all the other users do not transmit. In this case, we have 
L = A'bs and Hj 2i = H„. Moreover, all BS's and the nth user 
builds a single-user MIMO channel H„. Then, the maximum 
achievable rate is equal to the capacity for the MIMO channel 
H n with power Vl, which is given by (4) 

It£l= max {sTlog 



det (l + H„S^Ht)]} 



where det(-) generates the determinant of a matrix, Tr(-) eval- 
uates the trace of a matrix, and H^™- 1 is the covariance matrix 

(n) 

of s^' . Correspondingly, we get the maximum achievable 

L (n) 

effective capacity of the nth user, denoted by Cm ax , as follows: 

4"i=-ilog(E H {e-^}), (12) 

for n = 1,2,..., K mu . Furthermore, we define the effective- 
capacity fraction for the nth user as the ratio between the 
traffic loads and the maximum achievable effective capacity. 
Denoting the effective-capacity fraction by C„, we define 
C n = C„/Cmax- Note that C„ can be readily obtained off-line 
based on the statistical information of wireless channels, and 
thus can be used to design the BS selection algorithm during 
the data transmission process. For presentation convenience, 
we sort {C n }n=i in the decreasing order and denote the 
permuted version by {C^^j)}^^, where C n m > CV( 2 ) > 
• • • > CV(A' mu ) indicates the order from the higher priority to 
the lower priority. In the rest of this paper, we use the term 
of user Tiii) to denote the user associated with the ith largest 
effective-capacity fraction. 

Clearly, for a higher C n , the nth user needs more wireless 
resources to meet its QoS requirements. Thus, in order to 
satisfy the QoS requirements for all users, we assign higher 
BS-selection priority to the user with larger C n . Following 
this principle, we design the priority BS-selection algorithm 
to determine VLl in each fading state and provide the pseudo 
code in Fig. [2] In the pseudo code given by Fig. [2] we use 
temporary variables , J and ^ to denote the subsets of BS's 
which have been selected and which have not been selected, 
respectively. 

As shown in Fig. [2] in each fading state the BS-selection 
procedure starts with the selection for user 7r(l), who has the 
highest priority. After picking one BS for user 7r(l), we select 
one different BS for user ir(2). More generally, after selecting 
for user 7r(j), we choose one BS for user Tr(j + 1) from the 
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BS-subset which consists of the BS's that have not been 
selected. This procedure repeats until L BS's are selected. For 
user-7r(j)'s selection, we choose the BS with the maximum 
aggregate power gain over the subset "P, where 7 7r (j), m denotes 
the instantaneous aggregate power gain between user ir(j) and 
the mth BS (see Eq. (Q~|i for its definition). In addition, after 
user-7r(AT mu )'s selection, if the number of selected BS's is still 
smaller than L, we continue selecting one more BS for user 
7r(l), as shown in line 06 in Fig. [2] and repeat this iterative 
selection procedure until having selected L BS's. Clearly, users 
with higher priorities benefitted more from the above algorithm. 
Also note that the mobile users' priority order is determined 
by the effective-capacity fraction, which adapts to the mobile 
users' QoS requirements. 

A.2. The Principle of the Block Diagonalization Technique 

The block-diagonalization (BD) precoding techniques T\A\ 
have been widely used for MIMO transmissions because of 
its low complexity. In this section, we also apply the BD 
technique for our QoS-aware BS selection framework. For 
completeness of this paper, the principles of the BD technique 
are summarized as follows. 

Given transmission mode VJIl = (QlMl), the idea of block 
diagonalization lfl4l is to use a precoding matrix, denoted 
by T^, for the nth user's transmitted signal vector, where 
n = n u G Ul for some u, such that H^T^ = for all 
i satisfying i ^ n and i G Ul- By setting Sm^ = Smt , 
where s^ L is the nth user's data vector to be precoded by 
we can rewrite the received signal as 



(n) 



->(*) i An) _ f (n) , c (n) 



where r OTi = H^T^ . Under this strategy, the nth user's 
signal will not cause interferences to other active users. Ac- 
cordingly, the MIMO broadcast transmissions are virtually 
converted to U orthogonal MIMO channels with channel ma- 

trices { \ . Thus, the nth user's maximum achievable 

rate, denoted by (P^ 1 ), is equal to the capacity of the 

equivalent MIMO channel T m , as follows: 



max< BT log 



detfl + f^H^lr^ 



(n) 



(n)^ 



= V L for n G Ul, where H 



(13) 
is the 



subject to Tr(E 

covariance matrix of s^j and V^' denotes the power allocated 
for the nth user under mode DJIl- Correspondingly, we will set 



(n) 



the service rate R„ of the nth user equal to R^' {V L )• Note 

(n) L 

that may not exist, which then results in a service rate 

equal to 0. Also, we set R n = R$ l (V L n) ) = for n £ U L or 
L = 0. For the procedures to determine the precoding matrix 

(n) 

Tyji of the nth user, where n = n u G Ul for some u, please 
refer to lfl4l . 

A.3. Derivation of Active-User Set Ul 

Note that given Ql we ma y not be able to accommodate 
all users, because of the limited number transmit antennas. 



01. Let A := {1,2,. . . , K mu }, A : 

02. WHILE (A ^ 0) 
For all n 6 A 

Temporarily setting Ul 



and My, = 



fen 



Ml. 



03 
04 
05 



Get T 



(n) 



A U {«}. 



. Then, set 



InWrW 



Fixing T}£ 



J_T (n) 



conj 



-.(«) 
■m. 



or. 



where Tq ' 



06. 



07. 
08. 



09. 
10. 
LI. 



is 

given by Eq. (TTJ; (■ o ■) generates 
element-wise product between two matrices with the same size; 
conj(-) yields the element- wise conjugation. 
Set 



In 



if < VJ n < 
otherwise . 



i 



H 



-.(") 

■<m L 



END 

Select u such that for all j £ A, j / u, the following condition 
(72 > Ij) or (7a = 7j & user U has higher priority than user j) 
holds, where the priority order is determined in Section EV-Al l. 
IF ro 2 > 0, A := A U {«} and A := A - {«}; else BREAK. 
END 

Set U L := A. 



Fig. 3. Pseudo codes of the block-diagonalization based joint channel-priority 
algorithm to determine the active-user set Ul in each fading state. 



Although several algorithms for selecting active-user set have 
been proposed 11151 . |[16l . they cannot be applied in the frame- 
work of this paper, because the QoS provisioning for mobile 
users are not addressed those in these algorithms. Next, we 
determine Ul through a joint channel-priority method for active 
user selections. The pseudo code of this algorithm is provided 
in Fig. [3] 

In the joint channel-priority algorithm provided by Fig. [3] we 
iteratively select users one by one into the set Ul- In particular, 
we use variables A and A to represent the temporary sets of 
users which have and have not been selected, respectively. 
As shown in Fig. [3] lines 02 through 10 describe loops for 
iterative user selection, where we pick one user in each loop 
until all users are selected (i.e., A = 0) or no more user 
can be accommodated (examined by line 12). Within each 
loop, given the existing active-user set A we examine the 
channel quality of each user after BD. Specifically, we first 
get the BD precoding matrix of the nth user. Then, we derive 
w n , which is average channel-power-gain after BD over all 
transmit antennas, representing the average channel quality, and 
also obtain ||H["^r^^ HlV-Ms line 06, which characterizes the 
instantaneous channel quality. We further define a variable 7„, 
as shown in line 07, where 7„ = 1 and j n = indicate 
that the channel quality is above and below the average level, 
respectively. Obtaining j n , our selection criteria are as follows. 
First, we desire to select the user with higher j n , implying that 
this user's current channel is better compared with its statistical 
channel qualities, which will more efficiently use the system 
resources towards this user's QoS requirement. Second, if two 
users have the same j n , we will select the user with higher 
priority. Following this criterion, in line 08, we pick the unique 
user from A in the current loop, whose index is denoted u. 
However, if Wu = 0, implying the maximum achievable rate 
equal 0. As a result, no more user can be admitted, including 
the uth user. We will then terminate the loop, as shown in line 



TO APPEAR IN PROCEEDING OF IEEE INFOCOM 201 1 



6 



09, to finish the selection process. 



B. BS Selection in Single-User Transmission Modes 

For single-user transmission modes, we have Ul — {n}, 
n G {1, 2, . . . , K mu }. Thus, at any time instant, there is 
only one user receiving data from multiple BS's through a 
single-user MIMO channel Accordingly, the maximum 

(n) 

achievable rate for the nth user is equal to the capacity of _ff f2 

(n) 

with power Vl, which is denoted by -R^ . However, even for 
the single-user case, the complexity of high of choosing VIl 
to maximize the achievable data rate is too high, since we 
need to examine all ( K £ s ) combinations. Norm-based antenna 
selections have been demonstrated to be effective in achieving 
high system throughput with low complexity in centralized 
MIMO system 0, |6), which can be also extended to BS- 
selection in distributed MIMO system. Specifically, for the nth 
user with specified L in our framework, we select BS's with 
L largest aggregate channel power gain. As a result, in each 
fading state we have K^ s K mu single-user modes as candidates 
for distributed MIMO transmissions. Given the transmission 
mode with Q,l and the active user n, we will set the service 
rate R n equal to R$ . 



C. The Optimization Framework for Transmission Mode Selec- 
tion and Resource Allocation 



min (e h IV] L ( <p L + V? q L 
s.t: 1). <t>L + E E 9L,n = 1, VH 




L=0 L=l n=l 



2)- =Vl, VL, H; 

n=l 



(14) 
(15) 



,L=0 



EE 



1L.. 



< 0, Vn. 



(16) 



D. Derivations of the Optimal Solution of Problem Al 



DA. The Properties of R { $ L (V L n) ) 

Before solving Al, we first summarize the properties of 
Ryft CPl) determined by Eq. d!3t . Based on results in flU, 

--(*>) 

the nth user's MIMO channel 1^ (after BD) can be converted 
to z[ n ^ parallel Gaussian sub-channels, where Zj/ is the rank 



of r OT[ . Correspondingly, the zth sub-channel's SNR is equal 
to Ejf 1 ,, where the square root of Ejf 1 is the zth largest nonzero 



singular value of r^j' . The optimal power p^' z allocated to 
the zth sub-channel follows the water-filling allocation, which 



is equal to p£ z = — l/ e L™z] + ' w h ere [ - ] + — max{-,0} 



and fj^ is selected such that J2z=i P*lL = ^l' '■ Since T. 
has only Z^f 1 non-zero singular values, we define 1 / = oo 



for i = Z 



(n) 



We have derived candidate (SIlMl) m multi-user and 
single-user transmissions modes, respectively. We still need 
to answer how to allocate power over different mobile users 
in multi-user modes and which transmission mode will be 
eventually used for distributed MIMO transmissions. In this show that R%£ (V 
section, we employ the probabilistic transmission to determine 
finally selecting which transmission mode. Specifically, we 
use multi-user mode (Q,l,Ul) determined through algorithms 
given in Figs. [2] and [3] with a probability denoted by (j>L, 
L = 0, 1, 2, . . . , -Kbs; also, we use single-user mode with BS- 
subset cardinality L and Ul with a probability denoted by qL.n, 
L = 1, 2, . . . , K^, n = 1,2, ...,K mu . Note that O is the 
probability of the case that nothing is transmitted. Clearly, the 
sum over all qi,,n and 4>l must be equal to 1. For multi-user 
mode, we denote power allocated to the nth user in transmis- 
sion mode (Ql,Ul) by P L n \ while the total power constraint 
is given by Eq. (|3). For presentation convenience, we further 



(n)-i + 



(n) _ ,£,(71) 



1 and l/s£\ — for i = 0. We can further 



is a strictly concave function and 

(») (,,(»)) 



dV 



(n) 



BT 

„(") 



holds. Moreover, if /i L 
1,2, 



(n) 



,Zl',we get: 



1 / £ L,7' 1 / t LM 



(17) 



for i = 



>(n) _ 



E7 1 
1 = 1 _(") 



R$ L (V^) = BT log (nj =1 41) +S« log M £ 



,(n) 



(18) 
(19) 



define = (</>i, 2 , • 



fc.) and q = (q lf q 2 , . . . , q*r b J 
with q^ = (ql,i> Ql,2, ■ ■ ■ , q l , A^ mu ) to describe the probabilistic 
transmission policy; we also define V = yPi, "P 2 , ■ ■ ■ , 'PK bB ) 
with "Pl = (P^ j'Pi, , • • • , 7-£ ffmu - ) ) to characterize the power 
allocation policy in multi-user modes. Then, we formulate 
the following optimization problem Al to derive the efficient 
transmission-mode selection and the corresponding power al- 
location policy: 

Al: Joint BD-PT based BS-usage minimization 



D.2. The Optimal Solution to Al 

Theorem 1: The optimal solution for optimization problem 
Al, if existing, is given by 

1 + 



(n) J 



0, 



if n G Ul] 
if n ^ U L ; 



for all n, L, and H, where is the square of T^'s 
jth largest singular value; is the unique solution 



(20) 
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satisfying the following conditions: 

1 BT 

,,(") _ ( in.L \ WBTa n i* fM\ wn 

- \BT8 n \' n ) llj=l \ £ L,j) 

Vn, L, H. 



(») r- 



(21) 



The corresponding optimal PT policy is determined by 



11 



1 



(22) 



where is the indication function and tp* is defined as 



tp* = min <^ min {ip L } , min {^i,n} 



(23) 



with 



^ 0<£<X bs ; 
LG[l,JsT bs ], Vn; 



if multiple iPl's and/or ipL,n's all equal to tp*, the correspond- 
ing transmission modes will be allocated equal probability with 
the sum probability equal to 1. The variables {Xn}n=i ^ 
constants over all fading state; given {A* } n =i, Ch l ls selected 

to satisfy the equation J2n=i vl ')* = ^ L f° r a ^ ^ 
accordingly {A* }^=i need to be selected such that the equality 
of Eq. ( fl6l ) holds. 

Proof: We construct AVs Lagrangian function, denoted 
by Jai(4>, q, V; A, C H ). as 

l 7 j4 i(0 ) q,P;A,C H )=EH{JAi(0,q,P;A,CH)} (24) 

subject to X)l=o ^ + El=i En=i = !. where 



JAi(</),q,P;A,C H ) 






5^9i,n 1 - > >n / 

n=l 

L e "l.wiA i ;_|_ gi;ne n £ 



E E ^lj 



-e n c n 



(25) 



In Eqs. d24b-(|25l), A = (Ai, A 2 , . . . , A#- mu ) and A„'s are the 
Lagrangian multipliers associated with Eq. ( fT&t , which are 
constants overall fading states and satisfies A„ > 0; {Cu,l}l=i 
are the Lagrangian multipliers associated with Eq. ([TBI in each 
fading state, and Chx - (Ch,i, Ch,2, • • • , CH,K bs ). 

The optimization problem Al's Lagrangian dual func- 
tion |10|, denoted by 3ai(A, Ch)> is determined by 

5ai(A,Ch) - , min {^li(0,q,'P; A,C H )j 

= E H ( min {j A1 (0,q,P;A,C H )i}- (26) 

subject to J2l=o + Ei=i En=i = 1 for a11 H La- 
grangian duality theory shows [10] that 3Ui(A, Ch) is always 
a concave function, whose maximizer is upper-bounded by the 



optimum of Al. We then denote the maximizer of 3ai(A, Ch) 
by (A*,Ch)- Also, we denote by (0*, q* ,*P*) the minimizer 
to Eq. d26l >. which varies with (A, £). Then, given (A*,Ch)> 
we have 

(0*,q*) = argminiJ^q^A*,^)} 

arg min ^ ^ (j> L ip L + ^ E lL,n^L,n \ , (27) 
,L=o 1=1 n=l J 



for all H, where tp^ and V-L.n is defined in Theorem [TJ 
and equation (a) holds by applying Eq. d25l l and removing 
the terms independent of tp. Solving Eq. ( f27b subject to 
E£o^ + Ef=iE»=i u 9i,„ = 1, we obtain Eqs. 
d23l . If multiple ipi's and/or ipL,n$ all equal to i/j*, which 
happens with probability zero when H has continues CDF, 
how to allocate probabilities across these modes does not 
affect the eventual results. Therefore, without loss of generality 
we allocate the corresponding transmission modes with equal 
probability while keeping their sum equal to 1. 

It is clear that = for n 6 U L . Next, we 

consider n G Ul- Based on Eqs. (I22l-(l23l. the opportunity 
of transmitting the data in a fading state will be given to only 
one transmission mode. Moreover, given 4>l = 1 for some 
mode A4l, the power allocations for other mode do not affect 
the Lagrangian function. Therefore, "P* L needs to minimize 
Jai(4>, q, "P> A*, Ch) under <f>L — L 4>j — f° r all j ^ L, and 
QL.n = for all L,n. We denote J^i(0, q, T 7 ; A*, Ch) under 
this condition by Jai,l{"P> A*, Ch l)- Then, applying Eq. ( fTTb , 
taking the derivative of Jai,lCP> A*, Ch l) w i tn respect to 
(w.r.t.) V^, and letting the derivative equal to zero, we get 



> («) e -^<K n) ) =0 , Vn, 



L, H. (28) 



Deriving and applying Eq. dT8b . we obtain Eqs. (f20b- (f2TT) . 
We further define 



f n ((/>,q,T) =E H {Ef=o(^ 



+qL, n e 0nRtlL + T,j, j7 t„<lL,. 



which are the constraint functions on the left-hand sides of 
Eqs. ( TToT ) and $15[ , respectively. The Lagrangian duality prin- 
ciple [10 1 suggests that the optimal objective value L of Al 
satisfies: 



L >3ai(A*,Ch) 



(29) 



Also, fw., L {Vl) and / n (0*, q*, V*) are the subgra- 
dients JJO) of Zai{\ Ch) w - rt - Ch,l and A„, re- 
spectively. We can further prove that the subgradients 
f n {<j>*, q*,P*) and /h.l^I) of 3ai(A,Ch) ™y con- 
tinuously with (A,Ch)- Thus, Zai{\ Ch) is differentiable 
and we have 6Zai{\ Ch)/ 9a « = /«(</>*, Q*, ^*) and 
flXii(A, Ch)/5Chx = / H ,L(PI)5(H)dH, where 5 (H) is the 
probability density function (pdf) of H and dH denotes the 
integration variable. 

It is clear that if f n (<f>* ,q* ,V*) = (for all n) and 
fn,L(V* L ) = (for all L and H) hold, 3ai(A, Ch) attains its 
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maximum. Since /h,l('? > l) monotonically varies with Ch,l as 
observed from Eq. (l20b-(f2TT>. we can show that for any A, there 
exists a C h ,l resulting in /h,d(?1) = 0, L = 1, 2, . . . , K hs . 
This implies that Ch L must be selected such that the equal- 
ity holds in Eq. ( TT~5T > under A*. Due to the concavity of 
3ai(A, Ch)> <%(A, Ch)/^^k is a decreasing function of A„. 
Also, we can readily show that d$(\, C'u)/^^n\\ n =o > 0. 
Then, if there does not exist A such that <9J(A, C'n)/d^n 
for all n, we have A* — > oo for some nth user and 
d$(\, C,ii)ld\ n > always holds. For this case, we get 
L > 3(A*, Ch) ~^ °°> implying no feasible solution for Al. 

In contrast, if there exists A* such that &2ai (A* , Ch) / d\ n = 
for all n, we have Ch l = Ch l an d tne obtained (0*, q*, V*) 
is feasible to Al. Moreover, we get L* — 3ai(A*,Ch) w i m 
zero duality gap iflOl by examining Eq. (|24j, implying that 
(0*^*,^*) given by Eqs. (|20]i-(|23]i under A* and Ch is 
optimal solution of Al, and thus Theorem [TJ follows. ■ 

Note that there are no closed-form solutions for the optimal 
Lagrangian multipliers Ch ar, d ^* ■ 1° eacn fading state, Ch l 
needs to be selected to satisfy fn,h{'P* L ) = 0, as discussed in 
the proof of Theorem [TJ which can be conveniently determined 
through numerical searching method in that /m^CP^) varies 
monotonically with Ch.l- Moreover, we can determine Ch 
through maximizing the Lagrangian dual function 3(Chj^) 
by using the gradient descent algorithm. Due to the concavity 
of 3(Chi-^)' the gradient descent algorithm will converge 
with appropriately selected step size. If the gradient descent 
algorithm does not converge with A„ approaching infinity, the 
optimal solution does not exist for Al, as discussed in the 
proof of Theorem [TJ which implies that the current wireless 
resources cannot simultaneously support QoS requirements for 
all of current mobile users. 

E. Pure PT Based BS-Selection 

We further consider the BS-selection framework based on the 
PT-only approach for multiple access across mobile users. In 
this framework, the system only considers the Kb s K mu single- 
user modes derived in Section IIV-BI and the mode transmitting 
nothing as candidates for distributed MIMO transmissions. This 
PT-only based framework also uses probabilistic transmission 
to determine which transmission mode is used. Then, we can 
formulate the corresponding BS-usage minimization problem 
subject to the same power and QoS constraints as in prob- 
lem Al, where only the probability vector assigned for the 
-Kbs-Kmu + 1 candidate modes can be tuned to minimize the 
average BS-usage. The detailed problem descriptions and the 
corresponding optimal solution is omitted due to lack of space, 
but provided on-line in |18|. It is clear that this framework is 
easier to implement as compared to the joint BD-PT approach, 
but it can only support the lower traffic load. 

V. TDMA Based BS-Selection Scheme 

We next study the TDMA based BS-selection scheme. In 
the TDMA based BS-selection, we also apply the priority BS- 
selection algorithm given by Fig. [2] when the cardinality L of 
VLl is specified. Obtaining fi^, we further divide each time 



respectively. The nth user's time-slot length is set equal to T x 
tL,n for n — 1,2,..., K mu , where tL, n is the normalized time- 
slot length. Moreover, we still use the probabilistic transmission 
strategy across different VLl generated through Fig. [2] where 
the probability of using il^ to transmit data is equal to (p^. 
Then, we derive the TDMA based transmission policies through 
solving the following optimization problem Al. 
Al: TDMA based BS-usage minimization 



min {L\ — min < Eh 

(t,0) 1 ' 




s.t.: 1) 



E< 



VH, 



2). ^ = VH > £=1,2,.. 



(30) 
(31) 



3). E H 



< 0, Vn, (32) 



,L=0 



where 4> and t are functions of H. In particular, we have 

<t> = (00, 01, 02, • ■ • , 0K mu ), * — (*1,*2, • ■ ■ ,*Kbs)> and t L — 
(tL,l,tL,2, ■ ■ - ,tL,K b J- 

Theorem 2: Problem A2's optimal solution pair (t*, <fi*), if 
existing, is determined by 



1 



9 /? (n) 



log 



A n°n % 



H,L 



for all L, n, and H, and 



= J 1, if L = arg min \ i + ££™ X* n e 
0, otherwise, 



(33) 



(34) 



for all L and H, where <5Jj L under given {A* }^=i 15 deter- 
mined by satisfying 2~2 n =i 'in = 1» an d {^n}n=i nee ds to 
be selected such that the equality of Eq. ( [321 holds. 

Proof: The detailed proof of Theorem [2] is omitted due to 
lack of space, but is provided on-line in ifTHl . ■ 



VI. Simulation Evaluations 

We use simulations to evaluate the performances of our pro- 
posed QoS-aware BS selection schemes for distributed MIMO 
links. The BS's deployment and the mobile users' positions 
are shown in Fig. EJa), where A'bs = 5 and K mu — 3. We set 
T = 10 ms and B = 10 5 Hz. We further assume that all users 
have the same number of receive antennas, all distributed BS's 
have the same number of transmit antennas, and the incoming 
traffic loads for all users are equal. Furthermore, we employ 
the following average power propagation model. Specifically, 
the average received power gain h n _ m is equal to G/d^ TO , 
where d n , m is the distance between the nth mobile user and 
the mth BS, G is a constant factor, and n is the path loss 
exponent typically varying from 2 to 6 lfT2"l . Without loss of 
generality, we let "P rc f = 1 and select G such that h n m = dB 
a t d n ,m — 50 m. Also, we set a^ h = dB for evaluating of the 
average interfering range (see Section III-Q . 
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Fig. 4. (a) The deployment of BS's and the positions of mobile users, where -?ST mu = 3 and K^, B = 5. (b) Simulation results of the average BS usage L 
versus traffic load under the specified delay-QoS requirements, where £ n = 10~ 2 and = 500 ms for all n; Mm = 3; 

the average interfering range versus traffic load under the same system setup as in (b). 
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Fig. 5. (a) The average BS usage L versus traffic load under the specified delay-QoS requirements, where = 50 ms, £ n = 10~ 4 , and N n 

n; k = 1. (b) Average BS usage versus k, where M m = 5. (c) Average interfering range versus k, where the system setup is the same as in (b). 



2 for all 



Figures |4jb) and HJc) compare the average BS usage and 
interfering range as functions of the incoming traffic load 
among our derived QoS-aware BS-selection schemes, including 
the joint BD-PT, TDMA, and PT-only based schemes. Fig.gjb) 
shows that as the traffic load increases, all scheme's average BS 
usages become larger to satisfy the more stringent QoS require- 
ments. However, the TDMA and PT-only based schemes' BS 
usages increase much more rapidly than our proposed BD-PT 
based scheme. This is because block diagonalization for multi- 
user MIMO communications can effectively take advantage of 
space multiplexing in removing the cross-interferences among 
all mobile users, and thus can achieve high spectral efficiency 
and system throughput. We can further observe that when the 
traffic load gets lower (larger), the PT-only based scheme needs 
less (more) BS's to satisfy the specified QoS requirements, 
as compared to the TDMA based scheme. Fig. HJc) plots 
the average interfering range caused by distributed MIMO 
transmissions, which displays the similar results to Fig. Etb). 
This is expected because the total used power in each fading 
state linearly increases with the cardinality L of selected BS- 
subset, as shown in Section Hl-DI 

Figure |5ja) plots the average BS usage against traffic load 
with more stringent QoS constraints than the constraints used in 
Fig. IUb). Under these more stringent constraints, the PT-only 
based scheme cannot support the specified QoS requirements 
for the incoming traffics and thus are not plotted in Fig. 0a), 
which implies that the PT-only based scheme only works 



efficiently with loose QoS constraints. We can observe from 
Fig. [2a) that the BD-PT based scheme generally outperforms 
the TDMA based scheme in terms of requiring fewer BS's, 
especially when the traffic load is high. As shown in Fig. EJa), 
for traffic load higher than or equal to 500 Kbits/s, the BS 
usage of the TDMA based scheme will reach the upper-bound, 
which is equal to K^ s . This implies that all wireless resources 
have been used up while the specified QoS requirements for the 
incoming traffic still cannot be satisfied. In contrast, the BD- 
PT based scheme can clearly support even higher traffic load. 
An interesting observation is that the TDMA based scheme 
performs slightly better than the BD-PT based scheme, when 
the traffic load is low and the number of antennas per BS is 
small. This is because the advantage of BD technique can be 
effectively used when the spatial-multiplexing degree order is 
high. However, clearly the small number of transmit antennas 
can already successfully support smaller traffic load through 
TDMA strategy, while the BD in this case is not very effective 
due to the limited number of transmit antennas, implying 
insufficient freedom for spatial multiplexing. 

Figures |5jb) and |3c) depict the average BS usage and 
interfering range, respectively, versus the parameter n, where 
k is defined in Section III-D1 which is the power increasing 
rate with the number of BS's selected for distributed MIMO 
transmissions. We can see that the average BS usage and 
the interfering range of the BD-PT based scheme are much 
smaller than those of the TDMA based scheme. As shown in 
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Figs. [5jb) and |_Ic), the lower delay bound and the smaller 
violation probability threshold, implying more stringent delay- 
QoS requirements, cause more BS usage and thus larger inter- 
fering range. This is because in order to satisfy more stringent 
QoS requirements, more BS's need to get involved with the 
cooperative downlink transmissions to achieve the high system 
throughput for all mobile users. This also demonstrates that 
our proposed schemes can effectively adjust the transmission 
strategy to adapt to the specified QoS requirements. In addition, 
the average BS-usage is a decreasing function of n but the 
interfering range is an increasing function. This suggests that 
we can use more power to tradeoff the lower implementation 
complexity in distributed MIMO transmissions. 

VII. Conclusions 

We proposed the QoS-aware BS-selection schemes for the 
distributed wireless MIMO links, which aim at minimizing 
the BS usages and reducing the interfering range, while sat- 
isfying diverse statistical delay-QoS constraints over multiple 
mobile users. In particular, we developed the joint block- 
diagonalization and probabilistic-transmission based scheme, 
the TDMA based scheme, and the pure probabilistic- 
transmission based scheme, respectively, to implement efficient 
BS-selection and the corresponding resource allocation algo- 
rithms for QoS provisioning of mobile users. Simulation results 
show that the joint block-diagonalization and probabilistic- 
transmission based scheme generally outperforms the TDMA 
based and pure probabilistic-transmission based schemes in 
terms of requiring less BS's for data transmissions and de- 
creasing the interfering range caused to the entire wireless 
networks. Moreover, the TDMA and probabilistic-transmission 
based schemes is efficient when the traffic load is not heavy. 
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